Translational and rotational friction on a colloidal rod near a wall 
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We present particulate simulation results for translational and rotational friction components 
of a shish- kebab model of a colloidal rod with aspect ratio (length over diameter) L/D = 10 in 
the presence of a planar hard wall. Hydrodynamic interactions between rod and wall cause an 
overall enhancement of the friction tensor components. We find that the friction enhancements to 
reasonable approximation scale inversely linear with the closest distance d between the rod surface 
and the wall, for d in the range between D/8 and L. The dependence of the wall-induced friction 
on the angle between the long axis of the rod and the normal to the wall is studied and fitted with 
simple polynomials in cos 6. 
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I. INTRODUCTION 

A particle suspended in a fluid, moving in the vicin- 
ity of a stationary wall, feels a viscous drag force which 
is larger than the viscous drag it would experience in 
the bulk fluid [H-ll]- This may intuitively be understood 
by considering the special case of a particle moving to- 
wards (away from) a wall: fluid needs to be squeezed out 
(sucked into) the gap between the particle and the wall. 
Even when the particle is relatively far away from the wall 
the hindering effects of the wall are still felt through the 
long-ranged hydrodynamic inteactions. This has impor- 
tant consequences for practical applications where flow 
and time are issues. Especially for microfluidic applica- 
tions [3| , where large surface to volume ratios are encoun- 
tered, it is important to understand the fundamentals of 
near- wall dynamics. 

When dealing with colloidal particles random forces 
should also be taken into account [5| . The random forces 
are caused by temporary imbalances in the collisions with 
the solvent molecules, and lead to diffusive (Brownian) 
motion of the colloidal particles. The diffusive behaviour 
of nanometer to micrometer sized particles near walls is 
essential for the transient kinetics of phenomena such as 
wetting and particle deposition on a substrate 6] . 

The (anisotropic) diffusion tensor D of a colloidal par- 
ticle is related to the anisotropic friction tensor S by the 
generalised Einstein relation D = ksTS , where ksT is 
the thermal energy. The friction tensor in the presence of 
a stick boundary wall is difficult to obtain theoretically. 
Analytical expressions in the creeping flow limit (appli- 
cable to colloidal particles) are known, but are limited to 
the case of a spherical particle [2], [7[ or to a non-spherical 
particle whose major (hydrodynamic) axes are aligned 
with the wall and which is far removed from the wall 
0, Q. In the general case, particles are not aligned with 
the wall and/or may not be far removed from it. One 
then has to resort to experiment or numerical evaluation 
to obtain the friction or diffusion tensor. 

Experimentally, optical microsco py [9 |4l2l| ■ total inter- 
nal reflection microscopy (TIRM) [l3l fl4j , and evanes- 



cent wave dynamic light scattering (EWDLS) [y, [l5|, [l6[ 
have been used to determine the diffusivity of particles 
near a wall. The latter two experimental techniques use 
the short penetration depth of an evanescent wave under 
total internal reflection conditions, where in EWDLS this 
is combined with dynamic light scattering. In EWDLS 
the different components of the diffusion tensor may be 
obtained from the intensity time-autocorrelation, but 
this requires several careful theoretical interpretations 
16]. 

Numerical evaluation of the friction on a particle can 
be performed in several ways: by numerical summation of 
the forces due to a large number of Stokeslets distributed 
over the walls and surfaces of the particles 0, U3 , pos- 
sibly including image singularities to efficiently capture 
the effect of a planar wall [18l42C| , or by a multipole ex- 
pansion of the force densities induced on the spheres, also 
with an image representation to account for a planar wall 

In this paper we will present an alternative way to de- 
termine the friction on a colloidal particle, using molecu- 
lar dynamics simulations which explicitly include the sol- 
vent particles. Because of the large difference in length 
scales between a colloidal particle and a solvent molecule, 
it is impossible to perform such simulations in full atom- 
istic detail. Some form of coarse-graining is necessary. 
Here we choose the Stochastic Rotation Dynamics (SRD) 
method to effectively represent the solvent [23[ . The sol- 
vent interacts with walls and colloidal particles through 
excluded volume interactions [23 - |26| . Using this ap- 
proach, we determine the friction on a shish-kebab model 
of a rod of aspect ratio 10, i.e. 10 touching spheres on a 
straight line, as a function of distance to and angle (9) 
with a planar wall. We will show that the functional de- 
pendence of the wall-induced friction enhancement can 
be reasonably well described by the inverse of the clos- 
est distance d between rod surface and the wall (for d in 
the range between one eighth the rod diameter and the 
rod length), and an angular dependence which may be 
expressed as a simple polynomial in cos 9. These results 
serve as a first example to show how SRD simulations 
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may be used to determine the friction on particles of non- 
trivial shape in a non-trivial orientation with respect to 
confining boundaries. 

Stokesian dynamics codes, such as those developed by 
Brady and co-workers [20| and Cichocki and co-workers 
pll [22l | , have other sophisticated methods to determine 
the friction tensor. The results from such methods are 
generally more precise (if sufficient care is taken in its im- 
plementation and choice of parameters) than those from 
SRD simulations, because in SRD the solvent dynamics 
is stochastic and the resolution seems to be limited by the 
collision cell size. However, as we will show, in practice 
the resolution is better than a collision cell size, and the 
influence of stochasticy can be severely reduced by taking 
long time averages, leading to an acceptable precision for 
making predictions. The biggest advantage of SRD over 
Stokesian dynamics is the extreme simplicity of the im- 
plementation of the SRD method when non-trivial shapes 
and complex confining boundaries are involved. In Stoke- 
sian dynamics methods one has to deal with complicated 
multipole expansions and image representations. If the 
embedded particle shape is non-trivial, or if the confin- 
ing boundaries are not planar but more complex, they 
have to be represented by assemblies of different sized 
spheres. In contrast, in SRD, being essentially a molecu- 
lar dynamics technique with a simple additional rule for 
momentum exchange, all that is needed is a rule to de- 
termine when a solvent point particle overlaps with the 
embedded particle or wall. The SRD method then au- 
tomatically takes care of all hydrodynamic interactions 
between complex shapes. 

This paper is organised as follows. In section |TT] we 
give details of the simulations method and the constraint 
technique by which we determine the hydrodynamic fric- 
tion. In section Hill we validate the method by comparing 
simulations of a sphere near a planar wall with known 
analytical expressions. Then in section HVl we study the 
friction on a rod. In section [V] we conclude. 



II. METHOD 
A. Simulation details 

In Stochastic Rotation Dynamics (SRD) 23] a fluid is 
represented by Nf ideal particles of mass m. After prop- 
agating the particles for a time 5t c , the system is parti- 
tioned in cubic cells of volume <2q. The velocities relative 
to the centre-of-mass velocity of each separate cell are 
rotated over a fixed angle a around a random axis. This 
procedure conserves mass, momentum and energy, and 
yields the correct hydrodynamic (Navier-Stokes) equa- 
tions, including the effect of thermal noise 23]. The fluid 
particles only interact with each other through the rota- 
tion procedure, which can be viewed as a coarse-graining 
of particle collisions over time and space. 

To simulate the colloidal s phe res, we follow our earlier 
implementation described in [25| . Throughout this paper 



our results are described in units of SRD mass m, SRD 
cell size ao and thermal energy fc^T. The number density 
(average number of SRD particles per SRD cell) is fixed 
at 7 = 5, the rotation angle is a — n/2, and the collision 
interval St c = O.lto, with time units to = ^(m/fcsT) 1 / 2 ; 
this corresponds to a mean- free path of Af roo ~ O.lao. In 
our units these choices mean that the fluid viscosity takes 
the value 77 = 2.5m/ao£o an d the kinematic viscosity is 
v = 0.5ao/£o- The Schmidt number Sc, which measures 
the rate of momentum (vorticity) diffusion relative to the 
rate of mass transfer, is given by Sc = y/Df ~ 5, where 
Df is the fluid self-diffusion constant j.26-28]. In a gas 
Sc ~ 1, momentum is mainly transported by moving par- 
ticles, whereas in a liquid Sc is much larger and momen- 
tum is primarily transported by interparticle collisions. 
For our purposes it is only important that vorticity dif- 
fuses faster than the particles do. 

Stochastic stick boundaries are implemented as de- 
scribed in Ref. [25]. In short, SRD particles which over- 
lap with a wall or sphere are bounced back into the sol- 
vent with tangential and normal velocities from a thermal 
distribution. The change of momentum is used to calcu- 
late the force on the boundary. We note that despite the 
fact that the boundaries are taken into account through 
stochastic collision rules, the average effect is that of a 
classical stick boundary as often employed in (Stokesian) 
continuum mechanics. Because we will determine fric- 
tions by taking long time averages, the average flow ve- 
locities close to the boundaries will be effectively zero in 
all directions. In this work we set the sphere diameter to 
D = 8ao, which is sufficiently large to accurately resolve 
the hydrodynamic field to distances as small as D/16, as 
already shown in Refs. [25|, l26| . 

The method will be validated here again by compar- 
ing the friction between a sphere and a wall with known 
expressions from hydrodynamic theory. 

Walls are present at z — and z — L z , i.e. the wall 
normal n is in the z-direction. Simulations containing 
a single sphere are performed in a box of dimensions 
L x = L y — L z = 80ao, corresponding to 1QD in each 
direction. Simulations containing a rod with its longest 
axis along n are performed in a box of 10D x 10D x 20D. 
All other rod simulations are performed in a box of di- 
mensions 201? x 10D x 20D. The latter boxes contain 
approximately 10 7 SRD particles. Note that the finite 
box dimensions imply that there will still be an impor- 
tant self-interaction of the rod with its periodic images. 
Larger boxes are computationally too expensive. In this 
work we will assume that the dependence of the hydrody- 
namic friction on angle and distance to the wall is dom- 
inated by the presence of the wall itself and that self- 
interactions with periodic images lead to a simple overall 
multiplication factor of the friction. 

All simulations were run for a time of 5 x 10 5 io, i.e. 
5 million collision time steps. In CPU time this corre- 
sponds to about 4 weeks on a modern single core proces- 
sor. Such long run times are necessary to attain sufficient 
accuracy in the determination of the friction, the method 
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of which is explained in the following paragraph. 



B. Determination of the friction 

The translational friction tensor S transforms the 
translational velocity v of the centre-of-mass of an ob- 
ject to the friction force F it experiences. Similarly, the 
rotational friction tensor Z transforms the rotational ve- 
locity u> around the centre-of-mass to the friction torque 
T. In formula: 



F = -S-v 

T = -Z-lu. 



(1) 

(2) 



In our simulations, the translational friction tensor is de- 
termined without actually moving the rod by measuring 
the time correlation of the constraint force F c needed to 
keep the rod at a fixed position [2!| [3(| : 



1 



k B T 



dr{F c a {t Q + T)F c {t o )) 



(3) 



where a, ft £ {x, y, z} and the subscript to the pointy 
brackets indicates an average over many time origins 
to. Similarly, the rotational friction tensor is determined 
from the time correlation of the constraint torque T c 
needed to keep the rod at a fixed orientation: 



Z a /3 



1 



knT 



dr(T^t a + r)T c p (t )) t 



(4) 



It is important to realise that there are two sources of 
friction on a large particle [24| • The first comes from the 
hydrodynamic dissipation in the fluid induced by motion 
of the particle and may be obtained by solving Stokes' 
equation, integrating the solvent stress over the surface 
of the particle. For translation of a sphere of diameter D 
in an infinite fluid this yields the familiar (isotropic) fric- 
tion & = 37T7?J). The second contribution is the Enskog 
friction [3ll |32| , which is the friction that a large parti- 
cle would experience if it were dragged through a non- 
hydrodynamic ideal gas, i.e. through a gas where the ve- 
locities of the particles are uncorrelated in space and time 
and distributed according to the Maxwell-Boltzmann 
law. For a heavy sphere with diffusing boundaries (which 
randomly scatter colliding particles according to the 
Maxwell-Boltzmann law) the translational Enskog fric- 
tion is given by £ e = § ^2-nmk B T 1 D' 1 {\ + 2 X )/(1 + x), 
where y = 2/5 is the gyration ratio. It has been con- 
firmed [2^, [26|, [3II - I33I ] that the two sources of friction act 
in parallel, i.e. that the total friction £ is given by 



1 _ 1 1 



(5) 



For simplicity a scalar friction is shown here, but these 
results apply equally well to each component of the fric- 
tion tensor. The parallel addition may be rationalised as 



follows. When a large particle is forced to move through 
a sea of smaller particles, it can dissipate energy through 
two parallel channels: 1) by dragging itself through this 
sea of small particles, resulting in more large-small col- 
lisions, or 2) by setting up a flow field in the solvent, at 
the expense of viscous dissipation in the solvent but with 
the advantage that the sea of gas particles is co-moving 
near its surface, resulting in less large-small collisions. 

In a real (experimental) colloidal suspension both the 
solvent density and the range of the colloid-solvent inter- 
action are larger than simulated here, which is why at any 
given time many more solvent molecules are interacting 
with each colloid. This leads to an Enskog friction which 
is orders of magnitude larger than the hydrodynamic fric- 
tion (it is important to note that for a hard sphere in a 
point particle fluid the Enskog friction increases faster 
than the hydrodynamic friction with increasing sphere 
diameter, £ e oc D 2 and £h cx D). Hence for real meso- 
scopic particles the total friction £ is practically equal to 
the hydrodynamic friction In our simulations the En- 
skog friction on one sphere is only approximately 4 times 
larger than the hydrodynamic friction [261 ] . 

In our simulations, we can easily determine the En- 
skog friction. In many ways the hard-colloid and SRD 
fluid system is like a hard sphere system, in which case 
the general structural features of the force autocorrela- 
tion (or time dependent friction) are (a) an initial delta 
function contribution whose integral is the Enskog fric- 
tion and (b) a slower contribution associated with cor- 
related collisions and collective effects, which is neg- 
ative for low and intermediate particle densities [3l|. 
The Enskog friction can therefore be read off as the 
peak value in the running integral of Eq. ([3]), S a ^(i) = 

(1/fcnT) Jj dr (F£(to + r)F|(* )) , at very short times 

t. Note that for continuous interactions there is no clear- 
cut separation as described here [3l[ . Applying the above 
procedure to the case of a sphere yields an Enskog friction 
which is in good agreement with theoretical prediction, 
as we will show in section UTT1 Eq. is then inverted to 
Oi = !/(!/£ ~ V£e) to obtain the hydrodynamic friction 
coefficient which would be measured for a particle with 
orders of magnitude higher Enskog friction. 



C. Choice of system coordinates 

Looking at Fig. [T] it is clear that the cartesian coordi- 
nates x, y and z may not be the most natural coordinates 
to use when considering the friction on a rod. In this 
section we will introduce coordinates which are better 
adapted to the symmetry of the problem. We will show 
this for the case of the translational friction, but similar 
results will apply to the rotational friction. 

By symmetry, given the orientation of the rod in the 
xz-plane, the components 3^ and B, vz must be zero (and 
hence also S yx and 2 ZJ/ ). The remaining 4 independent 
components are the diagonals and the xz-component 
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FIG. 1: Rod and wall geometry. Z is the height of the 
rod centre-of-mass above the wall, d is the closest distance 
between the rod surface and the wall, which has its normal in 
the n direction. The unit vector ui is along the long axis of 
the rod, U2 is perpendicular to the rod and pointing parallel 
to the wall surface, and U3 is perpendicular to the previous 
two. The ange 8 is defined through n • ui = cos#. 



(which is equal to the zx-componcnt). Their relative im- 
portance in general depends on the orientation of the rod. 
In order to simplify this dependence as much as possible 
we transform the friction tensors to a rod-based orthog- 
onal coordinate system, defined as follows (see Fig. [IJ: 



ui 
u 2 

u 3 



u 

n x u 
|fl x u| 
u x (n x u) 
In x u| 



(6) 
(7) 

(8) 



In short, Ui is along the rod, u 2 is perpendicular to the 
rod but parallel to the wall, and U3 is perpendicular to 
the previous two. If 9 is the angle between the long axis 
of the rod and the wall normal (n • u = cos 9), with 
9 £ [0,7t/2], then the transformation matrix U between 
the cartesian frame and this new coordinate system is 
given by 



U 



sin 9 0— cos ( 

1 
cos 9 sin6> 



(9) 



With this we can calculate the transformed friction tensor 
as 



S = U T HU = 



cli 
e e± 2 



(10) 



where the friction components are given by 

£|| = sin 2 9 S xx + 2 sin 9 cos 9 E xz + cos 2 9 S zz (11) 
tu = 2 OT (12) 
C-L2 = cos 2 9 — 2 sin 9 cos 9 "i" sin 2 9 r E, zz (13) 
£' = sin0 cos (E zz - E xx ) + (sin 2 9 - cos 2 9) E xz 

(14) 

It will turn out that the mixing term £' is small relative 
to the other terms. 
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FIG. 2: (a) Autocorrelation of the constraint force needed to 
keep a sphere at fixed position Z — 3/4D above a wall. Forces 
perpendicular to the wall (zz) are shown in black; forces par- 
allel to the wall (xx) in red. (b) Running integral of the 
data in (a) (divided by ksT). The peak value at short times 
yields the Enskog friction (dashed line) £ e , which determines 
the short time self-diffusion of a colloid in a Brownian bath 
without hydrodynamic interactions. The total friction £ is es- 
timated from the long time limit of the running integral. The 
hydrodynamic friction £h is obtained by applying Eq. l[5]l. 



III. VALIDATION: FRICTION ON A SPHERE 
NEAR A WALL 

To validate our method, we first determine the friction 
on a single sphere as a function of its height z above a 
wall, and compare with known theoretical expressions. 

We have determined the constraint force autocorre- 
lations on a sphere of diameter D = 8ao for a series 
of heights ranging from z/D = 0.6 to 5.0 in a cubic 
box of size L — 10D along each axis. An example for 
z/D = 0.75 (i.e. with a gap width of 0.25D between the 
wall and the bottom of the sphere) is given in Fig. Uta) 
where we show both the perpendicular (zz) and parallel 
(xx) components. The structural features of the force au- 
tocorrelations are, as expected [3l| , an initial delta func- 
tion contribution, whose integral is the Enskog friction, 
and a slower contribution associated with correlated col- 
lisions and collective effects, which is negative for our 
relatively low particle density. 

The running integral of the constraint force autocor- 
relation (divided by fc^T) is shown in Fig. [2](b). The 
peak value at short times measures 0.74 x 10 3 in Simula- 
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FIG. 3: (a) Hydrodynamic friction enhancement due to the 
presence of a wall as a function of normalised height Z/D of 
the sphere's centre above the wall. Simulation results (circles 
and squares) are compared with theoretical expressions (black 
lines) for the friction for perpendicular motion towards and 
away from the wall, Eq. (|15p 2], and for motion parallel to the 
wall, Eq. (|16[) 1]. (b) Hydrodynamic friction enhancement 
as a function of normalised closest distance d/D, with d — 
Z — D/2, on a semi- logarithmic scale. The simulation results 
(circles and squares) are compared with lubrication theory 
asymptotes [2,13 which are valid for d <C D (red lines). 



tion units, which is in good agreement with the expected 
Enskog friction of £ e = 0.69 x 10 3 . We note that the 
measured peak value of 0.74 x 10 3 was found consistently 
for all distances between sphere and wall. This confirms 
that the value of the Enskog friction is a local effect, 
unaffected by the geometry of the surroundings of the 
sphere. 

After the Enskog peak, the running integrals converge 
slowly to their final values, in the particular case shown 
E zz = 0.35 x 10 3 and E xx = 0.23 x 10 3 , but with other 
values for other distances between sphere and wall. Using 
the measured Enskog friction, we then calculate the per- 
pendicular and parallel hydrodynamic frictions. For large 
values of z/D these hydrodynamic frictions converge to a 
value £°° = 230 ±20, which is in good agreement with the 
expected value of £ theor = ZirqD/{\ - 1A5D/L) = 220, 
where the factor between brackets takes into account fi- 
nite system size effects [26|, [U [35| . 

In Fig. |3{a) we plot the resulting perpendicular (cir- 
cles) and parallel (squares) hydrodynamic frictions, nor- 
malised by as a function of normalised height z/D 
of the sphere's centre above the wall. Clearly, the fric- 
tion is enhanced greatly as the sphere comes nearer to 
the wall. A theoretical derivation for the perpendicular 
friction enhancement A_l was given in 1961 by Brenner 
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(15) 



where a — cosh (2z/D). Fig. (3Ja) shows this theoret- 
ical result as a solid line. There is no exact analytical 
expression for the parallel friction enhancement Am. A 
commonly applied approximation due to Faxen [l| , which 
deviates less than 10% from the result of precise numer- 
ical calculations for z/D > 0.52 and gives essentially the 
same results for z/D > 0.7 0, is the following: 
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(16) 



Fig. Efa) shows this expression as a dashed line. 

When the closest distance d = Z— D/2 is much smaller 
than the sphere diameter, i.e. when d <C D, the Stokes 
equations can be solved asymptotically, leading to so- 
called lubrication forces 0, 0|- For the perpendicular 
friction enhancement the lubrication prediction diverges 
as the inverse closest distance 121: 



lim Aj_ = — . 

d/D->0 2d 



(17) 



For the parallel friction enhancement the lubrication pre- 
diction diverges logarithmically [3], 



lim An « 0.9588 



8 , 2d 
— m — , 
15 D 



(18) 



where the constant 0.9588 has been determined by fit- 
ting to precise numerical calculations [7] . An interesting 
question is whether our simulations are able to repro- 
duce these asymptotes. Lubrication forces are in essence 
a hydrodynamic effect and the SRD method in principle 
resolves fully the hydrodynamics, at least down to the 
scale of a cell size ao- The resolution is even better than 
this because of the averaging effect of the applied random 
grid shift [3a |. 

Fig. [3Jb) shows the friction enhancements measured in 
the simulations and the lubrication predictions against 
the logarithm of d/D, thus emphasising the behaviour at 
small distances. We observe that our method predicts 
a perpendicular friction which approaches the lubrica- 
tion prediction (solid red line) from above, but does not 
yet reach this limit within the range of distances stud- 
ied. The approach is quite slow, which is in agreement 
with the exact expression Eq. (fT5|) . The parallel friction 
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FIG. 4: Running integral of the autocorrelation of the con- 
straint force needed to keep a rod at fixed position Z — 4.75D 
and orientation 9 = 30°. The diagonal components represent 
the magnitude of the friction anti-parallel to the direction of 
motion, for motion along ui (£ii), U2 (£xi)> and U3 (£±2)- The 
mixing term £' represents friction along the ui direction for 
motion along the u.3 direction (and vice- versa). The mixing 
term is always much smaller than the diagonal components. 
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FIG. 5: Hydrodynamic translational friction components £|| 
(black circles), £xi (red squares) and £x2 (green diamonds) 
as a function of rod height Z, for four different angles 8. Solid 
lines are fits using Eq. (ffij)) . 



approaches the lubrication prediction (dashed red line) 
from above too, but here the lubrication limit is reached 
already for d/D < 0.1, in agreement with observations 
by Goldman et al. 0). 

The good agreement of our results with the theoreti- 
cal predictions for a sphere, both at small and large dis- 
tances, gives us confidence that the same method may 
also be applied to the case of a rod for which no theoreti- 
cal predictions are known. In the remainder of this paper 
we will focus on closest distances in the range of D/8 to 
10D. We will show that in this range of distances the 
wall-induced friction enhancement can be approximately 
described by additive contributions scaling linearly with 
D/d. We emphasise that this scaling is not exact, but 
serves to represent our measurements in a compact func- 
tional form which may be useful for future simulations. 
The similarity of the d _1 scaling with Eq. (IT71) is prob- 
ably coincidental because lubrication theory generally is 
not valid at such large distances. 



IV. FRICTION ON A ROD NEAR A WALL 
WITH L/D = 10 

A. Translational friction 

A typical example of the running integral of the con- 
straint force autocorrelation for a rod near a wall is given 
in Fig.fJ] The diagonal components £||, an d £x2 rep- 
resent the magnitude of the friction anti-parallel to the 
direction of motion, for motion along Ui, U2, and U3, re- 
spectively (note that in our simulations we do not really 



move the particles). The mixing term £' represents fric- 
tion along the Ui direction for motion along the U3 direc- 
tion (and vice- versa). The mixing term is always found 
to be at least one order of magnitude smaller than the 
three diagonal components, so to a first approximation 
may be neglected. The fact that the mixing term is al- 
ways much smaller than the diagonal components shows 
that iii , U2 and U3 are indeed close to the principal axes 
of the friction tensor. Moreover, the convergence of the 
friction data was confirmed by performing duplo runs for 
most systems, yielding identical results (including, for 
example, the oscillation visible in £xi(i) near t = 50to)- 
The Enskog friction was again determined from the 
peak value of the running integral at short times and the 
total friction from the limiting value at large times. The 
resulting hydrodynamic frictions are presented in Fig. [5] 
(symbols) as a function of distance between rod and wall 
for four values of the angle 6. We find that a reasonable 
approximation (within the range of distances studied) for 
the translational friction components is to treat the wall 
effect as additional to the bulk friction, with a dominant 
dependence on the inverse smallest distance between the 
surface of the rod and the wall and an angle-dependent 
prefactor. The smallest distance is defined as follows: 
a shish- kebab rod consisting of L/D spheres, with its 
centre-of-mass at height z — Z, making an angle 9 with 
the wall normal, will have a smallest distance d\ with the 
wall at z = given by 



d x =Z- 



L/D - 1 



cos f 



D. 



(19) 



Another wall is present at z 



with a smallest 
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FIG. 6: Prefactors A\\(0) (black circles), A±i(0) (red 
squares), and A±2(9) (green diamonds), as obtained from fits 
to translational friction coefficients using Eq. (|20p. versus rod 
angle 6 (inset) or cos 6 (main plot). Solid lines in the main 
plot are fits using Eqs. ([21]) - (f25)) . 



distance di to the surface of the rod given by d-i — 
L z + d\ — 22 . Within our approximation, the transla- 
tional friction components may be expressed as 



(20) 



where £ Q £j can be any of the £_u, or £j_2 compo- 
nents. The fits are presented in Fig. [5] as solid lines. 
The z — oo values and prefactors A are estimated by 
a least-squares-fit. Note that the prefactors A may in 
principle also depend on the aspect ratio p = L/D; 
our results apply to the case p = 10 only. We find 
the following results for the bulk values (in our units): 
£fj° = 860 ± 20, fft = 1250 ± 50 and fj° 2 = 1500 ± 80, 
independent of the particular value of 9. These values 
may be compared to the approximate theoretical predic- 
tions [37i £\\ eor = 2Trr]L/(\np - 0.207 + 0.908/p) = 575 

and £,f eor = Airr/L/ (lup + 0.839 + 0.185/p) = 800 valid 
for a cylindrical rod in an infinite solvent bath. The 
frictions we measure are higher because of unavoidable 
self-interactions between the rod and its periodic im- 
ages which overall tend to increase the friction. In other 
words, even when the walls are infinitely far apart (in the 
^-direction), the periodic boundaries in the other two di- 
rections still cause friction enhancements of each of the 
friction components. We find that the self-interactions 
enhance the component £fj° by the same amount, inde- 
pendent of the actual rod angle 9. Similarly the values for 
£li and £2° 2 are consistently the same for all rod angles. 

In Fig. [5] (inset) we present the prefactors A, as ob- 
tained from the fits, as a function of rod angle 9. In the 
absence of theoretical predictions we have tried several fit 
functions. Good single powerlaw fits can be made when 
the prefactors are plotted against cos 9, see Fig. [5] (main 
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FIG. 7: Hydrodynamic rotational friction components £|| 
(black circles), fxi (red squares) and fj_2 (green diamonds) 
as a function of rod height Z, for four different angles 8. Solid 
lines are fits using Eq. (^7)) . 



plot), resulting in the following fit functions: 

A\\(0) = A§ +B\\ (1/2 -cos*?) 2 (21) 
A X i(6) = A° L1 +B ±1 (l-cos9) 4 (22) 
A ±2 (0) = A° 2 +B ±2 (l-cos0) 4 , (23) 

with Am = 0.044±0.003, B n = 0.23±0.02, A° X1 = 0.064± 
0.005, B ±1 = 0.17 ±0.02, A° ±2 = 0.060 ± 0.004, and 
B ±2 = 0.31 ± 0.02. 

In summary, if only one wall is present near a rod of 
L/D =10 with closest distance d, the translational fric- 
tion components are approximately given by 



£11 

e± 2 



HI 



c « J 1 



0.044 + 0.23 (1/2 - cos 9) 2 ^ p ! ■ ) 



7> 



X 



0.064 + 0.17(1 - cos ( 



0.060 + 0.31 (l-cos( 



7> <*» 



B. Rotational friction 

Figure[7] (symbols) shows the hydrodynamic rotational 
friction coefficients as a function of distance between rod 
and wall for four different values of the angle 8. The 
component 0i (circles) represents the rotational friction 
for rotation around the long (iii) axis, (+1 represents the 
rotational friction for rotation around the U2 axis, and 
C_L2 the rotational friction for rotation around the U3 axis. 
The mixing term £' was again found to be at least one 
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order of magnitude smaller, and is therefore neglected in 
the following analysis. 

Similarly to the translational friction, we have treated 
the wall effect as additive to the bulk rotational friction, 
again with a dominant inverse dependence on the small- 
est distance d. Denoting rotational friction components 
with C, these may be expressed as 

l + C af} (9) ■ u D 



(a/3 ~ Cq 



di 



(27) 



where Cq/3 can be any of the £j_x, or £±2 components. 
The fits are represented in Fig. [7] as solid lines. Again 
note that the prefactors C may also depend on the aspect 
ratio p — L/D. We find the following values for the bulk 
values (in our units): Cm° = (0.22 ± 0.01) • 10 5 , Cxi = 
(7.4 ±0.1) • 10 5 and C12 = (6. 7 ±0.2) ■ 10 5 . Theoretically, 
the expression for a cylindrical rod of aspect ratio p = 10 
for the latter two reads [37j]: C± eor = ^L 3 /[i{\np - 
0.662 + 0.917/p)] = 7.7 -10 5 . We consider the rather good 
agreement with our results to be somewhat fortituous: 
the theoretical results have been derived for a cylinder of 
length L and diameter D in an infinite bath, whereas we 
have simulated a succession of spheres in a finite bath. 
Because forces on the the extremes of a rod have the most 
important contribution to the torque, the magnitude of 
the rotational friction on a rod is much more sensitive to 
the shape of its extremes than the translational friction. 
The rounded extremes of our model would correspond 
effectively to a cylinder of smaller length (for example 
for a cylinder with L = 9D a rotational friction of £j_ = 
6.0 • 10 5 would be predicted). 

The rotational friction around the long axis is less than 
3% of those around the two perpendicular axes, and rel- 
atively it remains much smaller also when the distance 
to a wall becomes very small. 

In Fig. [3] (inset) we present the prefactors C, as ob- 
tained from the fits, as a function of rod angle 6. Good 
single powerlaw fits can again be made with our measure- 
ments when they are plotted against cos# (main plot), 
resulting in the following fit functions: 

(28) 
(29) 
(30) 
1 — 



C\\(6) = C?,+E\\ (l-cos0) 4 



Cu(0) 
C ±2 (6) 



°_L1 
L- i 



cos 9) 4 



with C,9 = 0.036 ± 0.002, 



E±i(l 

,_1_ 2 -r E±2 (1 - COS0) 4 

£|| = 0.080 ± 0.006, 
0.094 ± 0.004, E±i = 0.74 ± 0.04, C\ 2 = 0.090 ± 0.004, 
and E ±2 = 0.097 ± 0.005. 

In summary, if one wall is present near a rod of 
L/D = 10 with closest distance d, the rotational friction 
components are approximately given by 

,\_ 1 • ■ - 

Ml 



Cn 



1 



0.036 ±0.08(1- cos 9f 



(31) 



Cn « Cx|l 

c± 2 « cr|i 



0.094 ±0.74(1 -cos 6>) 4 (32) 
0.090 ±0.097(1 - cos 9) 4 ^}(33) 
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FIG. 8: Prefactors C\\{9) (black circles), C ±1 (9) (red 
squares), and C±2{0) (green diamonds), as obtained from fits 
to rotational friction coefficients using Eq. Q27p . versus rod 
angle 9 (inset) or cos 9 (main plot). Solid lines in the main 
plot are fits using Eqs. p8| - (pO)l . 



V. CONCLUSION 

We have shown that Stochastic Rotation Dynamics 
simulations can be used to measure the hydrodynamic 
friction on an object by constraining its position and ori- 
entation and analysing the time correlation of the con- 
straint force. 

In this work we have applied the method to a rod of 
aspect ratio L/D = 10 near a wall. The main result is 
summarised in Eqs. $M$i - flM]) and d3U) - ([55]). Reason- 
ably good fits of both the translational and rotational 
friction could be made with the inverse of the closest dis- 
tance d between the rod and wall, at least in the range 
de[D/8,L]. 

We have found that for a noticeable friction increase 
the closest distance between rod and wall needs to be on 
the order of the rod diameter. Also, in agreement with 
common sense, the friction increase is strongest when the 
rod lies parallel to the wall (cos 9 — 0). For translations, 
the friction component £j_2 is the largest, as this corre- 
sponds (at least partially) to motion to and from the wall. 
For rotations, the friction component C±i is the largest, 
as this corresponds to motion to and from the wall of the 
extremities of the rod. 

In anticipation of an accurate theoretical treatment of 
this system, we have fitted the angular dependence of the 
friction components and found good fits in most cases 
with (1 — cos6>) 4 . We do not have a motivation for this 
functional form except that, intuitively, the friction in- 
crease must be relatively larger when a larger area of 
the rod is exposed close to the wall. Hence an increasing 
function of angle 9 is expected. The only exception seems 
to be the angular dependence of the parallel translational 
friction £| | , for which the smallest friction increase occurs 
at an intermediate angle of 60 degrees (see the inset of 
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Fig. in]). Because we cannot give a full physical motiva- 
tion, the scalings we have presented here are possibly not 
exact. However, they do serve to represent our measure- 
ments in a compact form which may be useful for future 
simulations. Such simulations are planned for the near 
future. 

More generally, the technique presented in this paper 
offers the posibility to determine with reasonable pre- 
cision the hydrodynamic frictions on complex objects, 
possibly with nearby complex boundaries, under circum- 
stances where theoretical calculations are too difficult to 



be performed. 
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